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Abstract 

We introduce an intrinsic notion of Holder-Zygmund regularity for 
Colombeau generalized functions. In case of embedded distributions be¬ 
longing to some Zygmund-Holder space this is shown to be consistent. The 
definition is motivated by the well-known use of Littlewood-Paley decom¬ 
position in characterizing Holder-Zygmund regularity for distributions. It 
is based on a simple interplay of differentiated convolution-mollification 
with wavelet transforms, which directly translates wavelet estimates into 
properties of the regularizations. Thus we obtain a scale of new sub¬ 
spaces of the Colombeau algebra. We investigate their basic properties 
and indicate first applications to differential equations whose coefficients 
are non-smooth but belong to some Holder-Zygmund class (distributional 
or generalized). In applications problems of this kind occur, for example, 
in seismology when Earth’s geological properties of fractal nature have 
to be taken into account while the initial data typically involve strong 
singularities. 

Keywords: Zygmund classes; Holder continuity; algebras of generalized func¬ 
tions; generalized solutions to differential equations 


1 Introduction 

When studying models of wave propagation in highly irregular media, e.g., in 
seismology, (hyperbolic) partial differential equations have to be considered with 
coefficients and initial data being generalized functions. The coefficients repre¬ 
sent the medium properties, which may be irregular, e.g., due to folds, fault 
zones, or junctions of different geological units as well as caused by long term 
physical processes within geological layers. Once the location of layer boundaries 
through jump discontinuities is completed, refined geological information is re¬ 
flected in a specific type of regularity patterns of the material properties within 
a certain unit. Often self-similar or multi-fractal behavior can be observed and 
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Holder continuity, and more generally, Holder-Zygmund spaces, were found to 
be a useful tool for a systematic qualitative analysis (cf. [HJ Chap.4] and 
Chap.IV] for a mathematical justification, and 00^1^11221 for seismological 
applications). 

In general, differential equations of the type mentioned above need not make 
sense or may fail to have solutions within the theory of distributions. However, 
embedding the singular coefficients first into an algebra of generalized functions, 
here Colombeau algebras, enables one to carry out a detailed analysis and yields 
unique solvability under mild conditions (cf. nnmaiini). 

A preliminary study of this procedure in Colombeau theory was undertaken 
in 0, where the focus was on microlocal properties and the regularization 
aspects of wavelet transforms. The feasibility of recovering Zygmund-Holder 
spaces of positive regularity in one space dimension after the embedding into 
Colombeau algebras was proven. In the present paper we extend this result 
to arbitrary dimension and regularity scale, although by slightly changing the 
definition proposed earlier. We also give first applications to simple differen¬ 
tial equations. In particular, we study a (l-l-l)-dimensional hyperbolic Cauchy 
problem with typical geophysical conditions on the coefficients. We show how 
the regularity of the measured wave depends on the regularity properties of the 
medium as well as of the initial value. 

The outline of this paper is as follows. After a brief introduction to the 
basics of Colombeau theory in Subsection 1.1 we devote Subsection 1.2 to a 
review of distributional Holder-Zygmund spaces and their characterization in 
terms of Littlewood-Paley decompositions and wavelet transforms. Section 2 
introduces the corresponding Colombeau-theoretic notion and discusses basic 
properties and illustrative examples. Section 3 presents simple case studies in 
applications to differential equations. 

1.1 Colombeau algebras of generalized functions 

We recall the basic facts about the so-called special Colombeau algebras on R". 
They can be defined on arbitrary open subsets, or even on smooth manifolds, 
contain the space of Schwartz distributions, and provide far reaching consis¬ 
tency with respect to analysis in distribution spaces. For further details and 
applications we refer to PEI El. 

The key ingredient of Colombeau algebras is regularization by nets of smooth 
functions and the use of asymptotic estimates with respect to the regularization 
parameter e. More precisely, it is based on a quotient construction as follows: 
we set (with I = (0,1]) 

Em := {(Me)5G/ G f I VAT CC R",Va G N” G N : 

sup \d°‘ue{x)\ = 0{e~^) as e ^ 0} 

xGK 

Af := {(Me)5G/ G £m I VAT CC R",Vm G N : 

sup |ue(a:)| = 0(e"^) as e ^ 0}. 

xGK 

Em is a differential algebras with component-wise operations, M is an ideal in 
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£m 5 and the special Colombeau algebra is defined as the quotient space 


g:=£y,/M. 


Since we consider only this type of algebras here we will omit the term ‘special’ 
henceforth. A representative of an element m of ^ will be denoted by (Me)e, 
and we will write u = [(ue)E] in this case. Smooth functions are embedded as a 
differential subalgebra simply by cr(/) = [(/)e]- 

To embed nonsmooth distributions we first have to fix a mollifier p € c5^(]R”) 
with unit integral satisfying the moment conditions J p{x) x°‘ dx = 0 when 
|a| > 1. Setting Pe{x) = e~^p{x/e), compactly supported distributions are 
embedded by io(w) = (w* pe)^ +Af. Using partitions of unity and suitable cut¬ 
off functions one may explicitly construct an embedding ip \ T)' ^ Q extending 
io, commuting with partial derivatives and its restriction to C°° agreeing with 
cr. Note that although ip depends on the choice of the mollifier p this rather 
reflects a fundamental property of nonlinear modeling where the interaction of 
singular objects depends on the regularization. Additional specifications of the 
regularization from a physical model may and should enter the mathematical 
theory at this point. 

The ring of generalized complex numbers C is defined as the set of moderate 
nets of numbers {{re)e G with |re| = 0{e~^) for some N) modulo negligible 
nets (|re| = 0(£™) for each m). 

1.2 Review: Holder-Zygmund regularity of temperate dis¬ 
tributions 

This section is a synthesis of related parts from the following sources: in the basic 
notation and setup of Zygmund spaces we stay close to [H]; all wavelet aspects 
are taken from [El; for further properties of Zygmund classes and related spaces 
we refer to EOIEI]. 

The result reviewed here is not new and neither are the techniques of its 
proof, given in the Appendix. However, we felt the need to unify various aspects 
which are crucial to our application later on. The concise summary of our efforts 
is the formulation of Theorem ^ 

Continuous Littlewood-Paley decomposition: Following [HI Sect.8.5] we 
introduce a continuous analog of the Littlewood-Paley decomposition. 

Choose If G real valued and symmetric such that |^| < 1 in supp((p) 

and (p{^) = 1 if |^| < 1/2. Put tp = |t=i= ’ grad(/3(^) so that the 

support of is contained in the annulus t/2 < |^| < t. Observe that we 

obtain a continuous partition of unity 



( 1 ) 


If / S is used as a Fourier multiplier for u G we will 

sometimes write this in pseudodifferential operator notation, i.e., f{D)u = 
* u where T and^ denote Fourier transform. 

Note that for any u G and T > 1 arbitrary we have 


(p{D)u + ip{D/t)udt/t = ip{D/T)u = T^{T ^ip){T.)*u 
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which converges to u in when T —> oo. This specifies the meaning of the 
following decomposition formula, which is valid in , 

/ OO 

•>p{D/t)udt/t. (2) 

Holder-Zygmund spaces: The classical Holder spaces for s > 0 

not integer, as well as their natural extension to s € N, the so-called Zygmund 
classes, appear in [HI Section 8.6] in an equivalent realization given by the spaces 
These are defined, for any real s, in terms of a continuous Littlewood- 
Paley decomposition by 

CKR”) := {u G =5^' I juj* := \\ipiD)u\\ sup {f\\i;{D/t)u\\< oo}. (3) 

t>i 

Let m G N. In the context of this paper we call a function g G o5^(R") a 
wavelet of (oscillation) order m if its first m moments vanish, that is 

J x°‘g{x) dx = 0 0 < jaj < m — 1 (4) 

and it is weakly radial 1 |15[ Chap. 1, Equ. (5.6)]), i.e., 

df 

X \9m\\ = i ve^o. (5) 

In particular, radial functions can always be normalized so that they satisfy ©■ 
We introduce the notation f{y) = f{—y) and fe{y) = s~^f{y/s) for a func¬ 
tion / on M" (and the bar denoting complex conjugation). If is a wavelet we 
consider the wavelet transform Wg : ^'(R”) —> C°°(R"' x R+), mapping u G =5^' 
into 

Wgu{x,e) = u * ge{x) V(a;, e) G R" x R+. (6) 

(Note that "denotes reflection, not inverse Fourier transform.) It is immediate 
that the image Wg{^') is contained in the subspace Om(R” x R+) of smooth 
functions all of whose derivatives have polynomial bounds in x, e and I/e (our 
notation deviates from (Sj where this space is denoted by o5^'(R" x R+).) On the 
space Om(R" x R_|_) we can define the wavelet synthesis operator Mg, mapping 
H G Om into an element MgH G oS^'(R”), defined by 

de 

MgH= lim / H{.,e)*g,— (7) 

r—*0M^oo S 

with convergence being understood weakly in o5^'(R") (cf. [Hi Chapter I, Sections 
24, 25, and 30]). With the aid of Mg distributions in can be reconstructed 
from their wavelet transforms modulo polynomials, i.e, for each u G there is 
a polynomial p on R" such that 

u = MgfWgu) +p. (8) 

The crucial observation that motivates the definition of Zygmund regularity 
within Colombeau generalized functions is a characterization which is valid for 
temperate distributions. As mentioned above this can be found in [H Chapter 
3] in the framework of Bony’s two-microlocal spaces. However we repeat the 
arguments given there in a ‘stripped down’ version appropriate for the current 
context. A detailed proof can be found in the Appendix. 
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Theorem 1. Let s be a real number and g G =5^(]R") have m vanishing moments. 
Let u be a temperate distribution on R”. 

(i) Let m > s. If m € (7|(]R") then its wavelet transform satisfies 

sup £"* ||ITgM(.,e)||^^ < oo. (9) 

ee(o,i] 

(ii) Let m > —s and g be weakly radial then 0 implies that there is uq G 

such that u - uo G 

Remark 2. Once more we want to emphasize that the statement of Theorem^] 
is included in the corresponding, and more general, results presented in Meyer’s 
book ca Chap. 3]. The characterization of Holder-Lipschitz-Zygmund regular¬ 
ity 0 < s < 1 via the asymptotic behavior of a wavelet-type transform at small 
scales has a forerunner in terms of Poisson integrals, e.g., in |‘24l VIL5] for the 
one-dimensional case and in Ha V.4.2] on M". 


2 Intrinsic Holder-Zygmund regularity of 
Colombeau functions 

2.1 Basic notions and coherence properties 

We recall that a mollifier is a function p G =5^(]R") with f p = 1. In addition, 
we will henceforth assume p to be radial. 

Mollifiers and wavelets: We restate the following facts from |a Sect. 3.3] 

(i) Let a G Nq with ja] > 1. Then the function := {d°^p)~ has m = jaj 
vanishing moments and for any u G c5^'(K”) 

d°'(u * Pe)(x) = £~^°‘^Wpc>u(x,£) V(x,£) G R" X R+. (10) 

In particular, (A^p)~ is a wavelet of oscillation order 2k. 

(ii) If / x°'p(x) dx = 0 when 1 < jaj < m — 1 then p. := —■^{pe) |e=i defines a 
wavelet of oscillation order m and for any u G 

dv 

u* pe{x) = u* p{x) + J Wpu{x,r) — V(a;,e) € M" X R+. (11) 

In view of Theorem^equation 11011 suggests to test for Zygmund regularity after 
embedding by looking at the asymptotic properties of high-order derivatives. 
The following definition is based on this idea and refines it in order to ensure 
mapping properties with respect to differentiations. Note that it differs from 
the definition proposed earlier in |^. 

Definition 3. Let s G R. and u = [(ue)^] G Q(R^). We say that u is of 
(generalized) Zygmund regularity s, denoted u G t/^(R"), if for a G Nq 

( 0(1) 0 < jaj < s 

l|a“welL=o = <^ 0(log(l/e)) |a|=seNo (£-0). (12) 

[o(£®“l“l) jaj > s 
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Remark 4. As a matter of fact, equation m and Theorem n directly suggest 
to include the third line in m of the above definition. This would be already 
suitable to characterize the embedded Zygmund classes (modulo smooth func¬ 
tions) among all embedded temperate distributions as can be seen from the 
proof of Theorem [3 below. However, if we want the family of spaces (s € K) 
to be a scale, in the sense that s' > s implies Q then the testing of 
decrease properties must not start at a derivative order which depends on the 
(prospective) regularity number. In particular, the case that s is an integer 
has to be taken into account, which is done here by the minimum possible, i.e., 
logarithmic, growth rate compatible with embeddings. 

Proposition 5. Let s, s', and r be real numbers. 

(i) The spaces are nested, i.e., s' > s implies C 

(ii) For each /3 G Nq we have a linear map ■. Ql ^ ■ 

(iii) Regularity of products: G^ ■ G^ ^ G* where p = r-|-sifr, s < 0, p = 
min(r, s) if max(r, s) > 0, and p = min(r, s)_ if max(r, s) = 0. (Here, c_ 
denotes any number c — a for tr > 0 arbitrary.) 

Proof. Part ('«_); If 0 < |q:| < s < s' the assertion is trivial. If |q:| = s < s' then 
||i9“w||/,cx, is 0(1) (s = s') or 0(log(l/e)) (s < s'), that is 0(log(l/e) in any 
case. 

The case |a| > s leaves us with three subcases for the asymptotic bounds 
of \\d°‘u\\^oc'. s < |a| < s' yields 0(1) which is 0(£®“l“l); s < |a| = s' gives 
0(log(l/e) and hence also 0(e®“l“l; finally, in case |a| > s' > s we obtain 
0(£®'“l“l being again 0(£®“l“L 

Part (ii): We use la with a replaced hy a + (3 and note that ja -I- /3| = 
|q;| -I- |/3|. This gives asymptotic bounds 0(1) if 0 < |q;| < s — |/3|, 0(log(l/£)) if 
joj = s — |/3|, and 0 (£®“l^l“|a| jf |q,| > § _ 

Part (iii): We may assume that r < s, the opposite case being completely 
analogous. Let u £ G(, v £ G),, and a G Nq. In estimating d°‘{uv) we use 
the Leibniz rule and thus have to find asymptotic upper bounds for the typical 
term of the form ■ d°‘~^Ve with (3 G Nq such that (3 < a. This is done by 
combination of the asymptotic growth information about each factor separately. 

If s < 0 then the largest growth is due to combinations of the form 0(£’’ 1^1) • 
Q(£s-|a|+|/3|) _ This proves the thirst case for the regularity p. 

If s = 0 we only have to check the case |q!| = 0 separately. To see this, note 
that adding —a in the exponents does not decrease the bounds established above 
and also captures any occurring logarithmic factors stemming from ||ue|j^oo. In 
order 0 the dominating terms are 0{£^) ■ 0(log(l/£)) = 0(£’’“°’) which proves 
the second case for p. 

Finally, let s > 0. Assuming |a| < s implies |/3| < r as well as |a —/3| < r < s 
and hence produces only 0(1) factors. If |a| = s then \a — j3\ = s if and only 
if \(3\ = 0 in which case the zero order bound for ||ue||^cx, is to be multiplied by 
log(l/£). Otherwise, i.e., if \(3\ > 0, then the factor corresponding to v gives 
only 0(1). It follows that we obtain the upper bound 0(log(l/£) if r = s > 0 
and if r < s. If |a| > s all possible nine combinations of 

upper bounds may have to be employed but 0{er — |q!|) is dominating all of 
them (since s > 0). □ 
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Remark 6. 

(i) Compare part (iii) of the Proposition with the distribution theoretic result 
on products in Zygmund spaces (cf. |H1 Prop.8.6.8]): If u G Cl, v £ Cl then u-v 
can be defined (as a weakly sequentially continuous bilinear map V xV ^ V) 
if r + s > 0 and gives an element of Zygmund regularity min(r, s). 

(ii) We note that the subalgebra Q°°, defined in [m Sect.25], reflects a some¬ 
what different concept of regularity. First of all, the ^““-property is tested on 
compact sets only with e-asymptotic constant with respect to derivative orders 
but dependent on the compact set. Furthermore, it is easy to give examples of 
Colombeau functions being very regular in one sense but not in the other: if p 
is a polynomial and x ^ smooth cutoff function then the class of x(a;)p(a:/e’') is 
in Q°° but it has poor Zygmund regularity if r > 0; on the other hand, for any 
s G M, e® sin(x/e) defines a -class which is not in Q^. 

Let p be a radial mollifier with all higher moments vanishing. (Hence p can 
be used to construct wavelets of any oscillation order.) Then we have the em¬ 
bedding 6p: 5^' ^ Q,v ^ [{v*Pe)e\- We show that under these embeddings the 
above definition of the subspaces I?* C ^ is compatible with the distributional 
Zygmund classes C®. 

Theorem 7. For any s G K.: 

(i) tp(C,^(IR")) C g-(R"). 

(ii) If u G ^'(R”) and tp(u) G ^J(R") then there is uq G C°°(R") such that 
u-uo G ^."(R”). 

Proof. Part (i): Let v & Cl and a G Nq. We work through all cases to be 
distinguished about the relation of jaj and s. 

laj > s and jal > 1: Application of ns and Theorem^] (i) (with m = jaj > s) 
yields |5“(w * p£)(a;)| = £:“l“l|Wp^w(x,e)| = (e —> 0) uniformly 

in X G M". 

0 < jal < s: In this case v G and we have ||9“(z; * Pe)\\L’=^ = || (9“u) * Ps\\l°° — 
||9“u||^oo||p|lii = 0(1), where we have used that O* C L°° if t > 0 1 ^11 
2.3.2/Rem.3]). 

jaj = s G Nq: Since d°‘v G 0° we obtain || Wp9“u(.,r)||= 0(1) in formula 
CD and hence 

f ^ fjj’ 1 

l|5“(u*p,)!!^,. < ||(a“u) *p||^„„ + c / - = o(iog(-)). 

jaj = 0 > s: Again by dJ and Theorem^] (i), noting that s < 0, we conclude 
that ^ 

\v * Pel < \v * p\ + C [ dr = 0(e®). 


Part (ii): Choose 2k > jsj. Then by m with p(^^) := (A^p)” and applying 
ca to Ve = V * pe we have 

The assertion follows from Theorem^ (ii) (with m = 2k > jsj). □ 
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The global L°°-bounds used in Definition 0 may be somewhat too restrictive 
in certain applications and instead of using a formulation like ‘is in modulo 
a very regular function’ we may prefer to use the following localized version of 
Zygmund regularity. 

Definition 8. Let il C R"- be open and s G R. The Colombeau function 
u = [(Me)e] G is Said to be locally of generalized Zygmund regularity s in 

D, denoted u G i^ all iG CC D and a G Nq 

f 0(1) 0 < |a| < s 

ll^“«ellL“(ic) =0(log(Ve)) |a| = sGNo (e - 0). (13) 

[o(e®“l‘^l) |a| > s 

2.2 Examples of regularity under composition 

Let (D) = risgRt/® [Qj.(n) denote the set of functions of arbitrarily high 
generalized local Zygmund regularity. In contrast to it we say that u has no 
Zygmund regularity, or regularity —oo, if it is not contained in UsgR^® (D). 

In the following we will consider the set Oc(K") of smooth functions all of 
whose derivatives are of the same polynomial growth, i.e., u G C°° and there is 
M G R such that for all a G Nq we have |u(a;)| = 0{\x\^) as \x\ oo; in this 
case, we will say that u is of growth order M. 

We determine the Zygmund regularity of a simple class of Colombeau func¬ 
tions obtained by scaling the arguments of smooth functions. These are not 
obtained by embedding of distributions and it is a special case of composing 
a smooth function with a generalized function. However, nontrivial regular¬ 
ity assertions about more general cases remain open at this stage. (A very 
useful result in Cf spaces concerning composition with smooth functions is |H1 
Prop.8.6.12].) 

Proposition 9. Let r be a real number. 

(i) Let / G Oc(R) of growth order M G R. Then Ui;{x) = fixje^) defines 
a Colombeau function u G t/(R) which is (at least) of local Zygmund 
regularity s if r < 1. We have s = —rM if 0 < r < 1 and M >0, s = l — r 
in case 0 < r < 1 and M < 0, and may put s = oo when r < 0. In general, 
we have no Zygmund regularity if r > 1. 

(ii) Let p be a polynomial of degree m 0. Then Ue{x) = p{xle^) defines a 
Colombeau function of local Zygmund regularity oo if r < 0. If r > 0 we 
have 

u G I/»joc(I^) ^ s < —rm. 


Proof. 

Part (i): If r > I we consider the (one dimensional) example Ue{x) = sin(a;/e'’). 
The derivative of order 2k, evaluated at a; = 7re’’/2, gives But this can 

never be dominated by for all fc G N and s fixed. Thus u has no Zygmund 

regularity. 

The other extreme case is r < 0 which always leads to £-independent bounds 
over compact sets in each derivative. Thus we have regularity of arbitrary order. 


We are left with the case 0 < r < 1. Let a G Ng then 

\d^u,ix)\ = £-’-i«i|(a“/)(^)l < + M)M_ 

Let X stay in a fixed compact set. If M > 0 the right-hand side is bounded by 
some constant times £:-’'(l“l+^) = Finally, if M < 0 all we can 

say (in general) is that the right-hand side is 0(£“’'l“l) = which 

is 0(1) if |q;| = 0 and 0(£^“’’“l“l) otherwise. 

Part (ii): The case r < 0 is obvious since all derivatives have upper bounds 
independent of e then. So we assume r > 0 and note that p is not the zero 
polynomial since it has degree m>\. 

Let a e Nq and assume 0 < |q;| < m, all higher derivatives vanish. We have 
\d°"Ue{x)\ = £“’’l“l|(i9“p)(a;/£’’)| which is e“’’l“IC)(e“’’(™“l“l)) = 0(e~™) if x 
varies in a compact set. Furthermore, since d°‘p is a polynomial (nonzero for 
some a of each occurring order) the estimates cannot be improved. 

Assume that u G Ql Since —rm is strictly negative is never dom¬ 
inated by a constant or logarithmic growth. Hence we have the conditions 
s — k < —rm when 0 < k < m. Setting k = 0 yields s < —rm. 

On the other hand, s < —rm is sufficient to establish the corresponding 
Zygmund regularity by the above estimates. □ 

We end this section with two examples falling into the range of the above 
proposition and that further illustrate the different behavior of the notions 
of Zygmund- and ^““-regularity, in particular, with respect to stability under 
smooth compositions. 

Example 10. 

(i) We have v = [{x/e)e] G and v S s < —1 (put m = r = 1 in the 

proposition above). Consider the composition u = sinov then u ^ Q°° but 
ugG^ iM = 0,r = 1). 

(ii) Similarly, v = [(1 -f x^ je)^ G G°° and r; G s < —1 (use the proposition 

with m = 2, r = 1/2). Since Vg > 1 for all £ > 0 we may form u = Ijv G Q. 
We observe that u ^ at x = 0, the values of the derivatives can be read 
off the coefficients in the power series expansion valid in the 

interval (—-^/e). From the proposition, with M = —2, r = 1/2, we deduce 
that u G . 


3 Application to linear differential equations with 
nonsmooth coefficients 

3.1 Solutions with classical Holder continuity 

We start with the simplest possible case of a differential equation and men¬ 
tion the well-known elliptic case only briefly. Finally, we sketch how a gain of 
regularity can be observed in the hyperbolic case too. 


9 



Primitive functions in one dimension: Let s be any real number and 
u G If u G I?'(R) is a primitive distribution of u, i.e., v' = u, then there 

is / G C°“(]R) such that 

v-feCt+^{R). (14) 

To see this we can employ an explicit parametrix of given as pseudodiffer¬ 
ential operator with symbol h{^) = where x ^ C“(R) vanishes near 

^ = 0 and xiO = 1 when |^| > 1. (Note that u G =5^'(R) and T{{h{D)uy — u) = 
(x — l)w has compact support; hence {h{D)u)' — m is smooth). It follows that 
{v — h{D)uy = u — {h{D)uy is smooth and so v — h[D)u must be. But h{D) 
being of order —1 maps C® into (7®+^ (see [HI Thm.8.6.14]) which proves (tCTl . 
Alternatively, we could state that v is locally in C®’*'^ in the sense that Lpv 
belongs to this space for any test function ip G T). 

Elliptic partial differential operators: Consider P{x, D)u = / where P 
is an elliptic partial differential operator of order m with coefficients and right- 
hand side / in C®, s > 0. Then u G (7®+™, i.e., we observe a gain in regularity 
by the order of the operator. More precise statements and related results can 
be found in Ch.3], a concise summary is |7| Thm.lT.l.T]. 

The embryonic hyperbolic case: As a resemblance of more realistic models 
from geophysics we consider the Cauchy problem 

dtU + a{x)dxU = 0,u\t=o=b (15) 

where a G (7®)®), 0 < s < 1, and b G (7®+^(M). In addition, we make the 
following strong positivity and boundedness assumption on the coefficient: there 
exist constants ci, C 2 such that 

0 < Cl < a{x) < C 2 for all a: G R. (16) 

This condition is justified, e.g., if a is of the nature of sound speed in a certain 
medium or fluid. 

The Cauchy problem US) is easily solved by the method of characteristics. 
We point out that, by continuity and positivity of the coefficient a, the char¬ 
acteristic ODE has indeed a unique C^ solution. To make this more explicit 
we define A{x) = dr/a(r). Note that A is C^, strictly monotone, and that 
|A(a::)| < |a:|/ci. Then we set 

u(x,t) = b(A~^(A(x) — t)) (17) 

which is directly checked to be the C^ solution of (CHI). As an introduction 
to the subject of the following two sections we investigate its Holder-Zygmund 
regularity in some detail. 

Proposition 11. Let u be the solution of 03) given by 03. Then the first 
order derivatives of u are Holder continuous of order s. 

Proof. Note that 1/a is in (7® which can be seen directly or, alternatively, be 
deduced from (HI Prop.8.6.12] since a is bounded away from zero. We proceed 
straightforward in two steps. 
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The function h{x,t) = A~^{A{x) — t) clearly is C^. We first show that 
its first order derivatives are Holder continuous with exponent s. We have 
gradh(x,t) = a(h(x,t)) ■ (l/a(a:), —1), which is bounded, and 

|grad/i(a:, t) — gradh{y, r)| 

< |a(h(a;,t))||(-^ - -^,0)| + |a(/i(a;, t)) - a(h(y,r))|-1)| 
a{x} a(y) a[y) 

< C{\x - y\" + \h{x,t) - h{y,r)\‘') < C'{\x - j/l® + |(a; - y,t-r)\‘‘) 

with generic constants C, C depending on a only. Hence grad/i is Holder 
continuous of order s. 

The second step is the composition with b. We have gradu = b'{h) ■ grad/i 
and therefore obtain 

|gradM(a;, t) — gradu(j/, r)| 

< |6'(/i(a;,t))||grad/i(a;,t) - gradI + \b'{h{x,t)) - b'{h{y,r))\\gradh{y,r)\ 

< C{\x - J/I® + \h{x,t) - < C'{\x - 2/1® + |(a; - 2 /,t- r)|®) 

where we have used the Holder continuity, as well as the boundedness, of grad h 
and h'. □ 


3.2 Primitive functions and a linear first order ODE 


The simplest inhomogeneous case is that of primitive functions in one dimension. 
Unlike primitive distributions, a Colombeau primitive function need not gain 
regularity, as the following examples illustrate. 

Example 12. The generalized constants [(l/e’')^], r > 0, do not have Zygmund 
regularity higher than —r but nevertheless are primitive functions of 0. As a 
consequence, any Colombeau function allows for primitive functions with Zyg¬ 
mund regularity arbitrarily low. Furthermore, all primitive functions of [(l/e’')£] 
are of the form [{x/ -I- c where c is any generalized constant. The latter can 
never be of Zygmund regularity higher than —r thereby showing the existence of 
Colombeau functions possessing no primitive function of any higher regularity. 


However, saving a minimum of the classical intuition, we can still control 
the regularity of primitive functions obtained from embedded distributions via 
integration. 

Proposition 13. Let u G 6p(C'®(M)), a;o G R arbitrary, and define a primitive 
function v by the representative Vg(x) = Ue{y) dy. Then v belongs to 

Proof. There is uq G Cf such that u = tp{uo). By itTHl we can find w G of 

uq such that we have ip{uo) = ipiw') + a{g) for some smooth function g. Hence 
there is {ne)e G Af such that 


(x) = / w'ps{y) dy + / f{y)dy + ns{x). 


We observe that, in general, any derivative of order I > 1 has the asserted 
asymptotic estimates since vi^\x) = ui so only the zero order estimate has 
to be investigated separately. 
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Using w' * Pe{y) dy = w * ps{x) — w * Pe{xo) we obtain, for any compact 
interval / containing x, xq, and of length |/|, 


|ne(a:)| < 2sup |w * pe(?/)| + sup (|/||/(y)| + \ne{y)\). 

yei yei 


The second term on the right-hand side is 0(1) on compact subsets with respect 
to X. Finally, since w G we deduce from Theorem |3 the required growth 

properties, according to regularity s -f 1, of the complete expression. □ 


In the proposition to follow we give a lower bound for the regularity of 
the solution to a linear homogeneous ODE with coefficient from a generalized 
Zygmund class. We will impose an additional boundedness condition on this 
coefficient and recall: v G G is said to be locally bounded if VK CC K" there is 
C, So > 0 such that sup,j,g^ < U for all 0 < e < Eq- 

Proposition 14. Assume s > — 1 and let a G G^ ioc(®^) such that Re(a) is locally 
bounded. Let 6 be a generalized constant which, considered as a generalized 
function, is of generalized Zygmund regularity t (t G K.). Then the unique 
solution u G ^(K) to the initial value problem 



a{x)u(x), m(0) = b 


(18) 


belongs to Gl^iod^) with r = s-|-lift>0 and r = t if t < 0. When t = 0 we 
have r = 0_ if s = — 1 and r = 0ifs>—1. Here, 0_ stands for any negative 
number, arbitrarily close to 0. 

Proof. Existence and uniqueness of the solution u follows from . A represen¬ 
tative is given by 

udx) =bee^o‘^-(y)dy 

where {bg)e is a representative of b. By our assumption on a we have on any 
compact set K 

= 0{\be\) (£^ 0 ). 

To find sharp asymptotic bounds for the derivatives we first investigate their 
algebraic structure. The following assertion is easily proved using the ODE 

(k) 

itself and induction on the derivative order k. Let fc G N then is a linear 
combination of terms of the following form: with to S N, 1 < m < /c, and 
A G N™ such that |A| = fc — to we have the expression 

m 

(19) 

f=i 

As noted above the first factor, u^, is Od^el), so we focus on the product of 
derivatives of ac¬ 
claim: for any s > — 1 we have, with the notation as in m, 

m f 0(1) fc < S -|- 1 

nil«^'^llL“(ic) = <^0(log(l/£)) k = s + l (20) 

[o(e*+i-'=) fc>s-hl 


on compact sets with respect to x. 
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• If s < 0 then each Aj > 0 > s and hence we have the asymptotic bound 
0(e’^«-l^l) = o(e™(®+i)-fc) = 0(e«+i-'=). 

• If s = 0 let n be the number of j’s such that Xj = 0. Then we have the 

asymptotic upper bound involving (log(I/e))”e“l''‘l = (log(I/e))”e'"“^. When 
m > 2 the second factor is where £ can compensate for the logarithmic 

terms. Hence we have a bound 0(£^“^). When m = 1 we obtain 0(log(l/£)) if 
k = 1 and 0(£^“^) otherwise (since n = 0 then). 

• Finally, we have to consider the case s > 0. We have to further distinguish 
three subcases for the relation between k and s + 1. 

Subcase fc < s + I: Since |A|=fc — m<A: — l<swe have that each Xj < s 
and hence an upper bound 0(1). 

Subcase A: = s + 1: Now |A| < s and for at most one j we have Xj = s, all 
others are less than s; hence we obtain an estimate of the form 0(log(l/£). 

Subcase/c > s + 1: Denote by n the number of j’s such that Xj = s and define 
N' := {j I Xj > s}, n' := |IV'| (the cardinality of N'). Put A' = 0 if j ^ W and 
X'j = Xj otherwise. The asymptotic upper bound in question is now expressible 
as 0((log(l/£))"£"'®“l'^'l). If n' = |A'| = 0 this clearly is 0(e®+^“^), so we may 
assume that n' > 1. Since fc — to = |A| > ns + |A'| we obtain £“1'*' I < 

Inserting this into the expression for the asymptotic upper bound we arrive at 
0((log(l/£)"£"'’)0(£” *+"*-'=). Here, the first factor is 0(1) since s > 0 and the 
second factor is 0(£®'''^“^), due to n' > 1 and to > 1, as claimed. 

Now we come back to m and use the information from 12011 . If t < 0 the 
order zero estimate implies r < t. By combining m with we see that 
regularity r = t can indeed be established for any s > — 1. 

If t = 0 the order zero estimate is logarithmic, due to |6e|, and forces r < 0. 
If s + 1 > 0 it is seen from (EDI) that r = 0 holds. In case s = — 1 and fc > 0 we 
have to cope with appearing upper bounds of the form 0(log(l/£)£“^). This 
requires subtraction of an arbitrary small, but still positive, number a in the 
exponent to incorporate the additional logarithmic factor. (Compare with the 
situation in the general multiplication result.) 

Finally, if t > 0 the factor |6e| = 0(1) and the regularity r = s + 1 is 
established directly from EOJ. □ 

Remark 15. 

(i) Note that if k is very large in it may happen that each Xj > s. In this 
case, a general upper bound will be of the form 0(£™®“l'^l) = 0(£™^®'''^^“*^). 
Since to may also become arbitrarily large this indicates that the condition 
s + 1 > 0 cannot be dropped in general while expecting Zygmund regularity 
of the solution. This is illustrated by the constant coefficient problem with 
ae{x) = ije^, 6=1 and r > 0. The solution (representative) is then exp(ja;/£’'), 
a sort of ‘standard counter example’ in Colombeau regularity theory. 

(ii) The boundedness condition on the real part of the coefficient cannot be 
dropped. Indeed, this can be seen from the constant coefficient problem with 
ae{x) = log(l/£) and 6=1. A Colombeau solution representative is given 
by exp(a;log(l/£)) which is not Zygmund-regular: The L°°-norm taken over a 
compact set K grows like £“"*('^1 if m{K) denotes the maximum of K. 
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3.3 A linear hyperbolic Cauchy problem 

As we have indicated in the introduction, if we think of modeling seismic wave 
propagation we may encounter fractal-like variations in sound speed, for ex¬ 
ample. By the very nature of coefficients representing physical observables like 
sound speed, density, elasticity tensors, we see that a positivity condition on the 
coefficient(s) is not artificial. We state a first regularity result for a simple model 
of this type in one space dimension. It hts nicely with the classical embryonic 
case discussed in Subsection 4.1. 

Theorem 16. Let a = [(uele] G ioc(®)i s > 0, and assume there are positive 
constants Ci, C 2 such that ci < ae{x) < C 2 for all a; G R and £ G (0,1). Let t 
be a real number and b £ Ql If u is the (unique) solution of the Cauchy 

problem 

dtu + a{x)dxU = 0, u{0) = b (21) 

then u £ Ql ioc(®^) with r = min(t, 1)_ if s = 0, and r = min(t, s -I- 1) if s > 0. 
(As above, min(t, 1)_ denotes any number approximating min(t, 1) from below.) 

Proof. We have to determine asymptotic upper bounds of all derivatives of 
Ue{x,t) = b^{A~^{As{x) — t)) on compact sets. 

We first note that the assumptions on a imply that he{x,t) = — 

maps a compact subset K of R.^ into a fixed compact subset K' of R, indepen¬ 
dently of e. Therefore when doing estimates on K we may essentially ignore the 
argument hg(xA) whenever appearing as inner function in compositions and 
write instead the supremum over K'. However, the chain rule will bring out 
derivatives of hg as additional factors. 

Thus the order zero estimate for is simply 

l|w£|lL”(if) - (22) 

In the following, let a G Ng such that |a| > 1. 

As a preparation we have to investigate the structure of the higher order 
derivatives of u^. To simplify notation we drop the subscript e in doing this 
algebra. 

Claim 1: d°‘u is a linear combination of terms of the following form: 

m n 

b‘^^\h{x,t)) ■ f)) ' {x) (23) 

i=i j=i 

where 1 < I < m < |Q;|,0<n< |q;|,0<A:< |q:|, |A| = m —I, and |/r| = |a| — to, 
with the notation A := (Ai)(T^ and /r := (/r^-)”^^. 

We prove (EH by induction on |a|. Concerning the inner derivatives when 
applying the chain rule we note that, by definition of A, we have dth{xA) = 
—a{h{x,t)) and dxh{x,t) = a(h{x,t))/o.{x). 

The base cases correspond to the first order derivatives dtu{x, t) = —b'{h{x, t))- 
a(h{x,t)) and dxu{x,t) = bfh{x,t))a(h(xA))/a(x)^ both complying with the 
structure of E3- 

Assume the claim to be proven already for |a| and let /3 G Ng with |/3| = 
|a| -I- 1. We distinguish the two cases (} = a + ei and (} = a + 02 {oj denoting 
the standard unit vector in direction j). 
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Case (3 = a + ei: By the induction hypothesis, d^u = dx{d°'u) is a linear 
combination of terms 

m n 

dx(b^'^Hh{x,t)) • /a!"{x)^ . 

i=i j=i 

Application of the Leibniz and chain rules yields four types of terms. 

Type 1 is 

m n 

b^^'^^\h{x, t j) a{h{x, t j) a{x)~^ ■ a^^'\h{x,t)) ■ a^^^\x)/a!^{x) 

i=i j=i 

which matches the claim with new quantities / + !, m+l, k + 1, and Am+i := 0 
Type 2, for any 1 < r < m, is 

and satisfies ll^ with fc + 1, m + 1, + 1, and Am+i := 0 instead. 

Type 3, for any 1 < r < n, is 

m n 

i=l j/r,j = l 

where we may use the new component + 1 in (P| . 

Type 4 is 



it^j)(^x) 

a^{x) 


i=i j=i 

and matches the claim with new quantities fc + 1, n + 1, and /in+i := 1. 

Case (3 = a + e 2 '- By the induction hypothesis, d^u = dt{d°‘u) is a linear 
combination of terms 

m n 

dt{b^^\h{x,t)) • ■X{a^^^\x)/a\x)). 

i=i j=i 

Application of the Leibniz and chain rules yields two types of terms. 

Type 1 is 


— a{h{x,t)) ■ a^^^\x)/{x) 

i=l j=l 

which matches the claim with new quantities 1 + 1, m + 1, and Am+i := 0 
Type 2 is 

m n 

—b^^\h{x, t)) ■ a^^^\h{x, t)) ■ a^^''^^\h{x,t)) a{h{x,t)) ■ a^^^\x)/{x) 

and satisfies (1^ with new values m + 1, A^ + 1, and Am+i := 0. 

The claim is proved. 
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According to claim 1 and the remark at the beginning of this proof we deduce 
that on any compact set we have 

m n 

Wd^u^h-iK) = o{\m\L-iKn ■ n ■ n (24) 

i=l i=l 

To evaluate this carefully we first focus on all the factors having bounds de¬ 
pending on a or its derivatives. 

With the notation of GSl define the sets Lq = {i \ Xi = s}, Li = {i \ Xi > s}, 

= {j I Mi = s}, and Mi = {j | Xj > s}. Let Iq = \Lo\ and define similarly 
h, mo, mi as the respective cardinalities. 

Claim 2: On compact sets we can give asymptotic upper bounds of the 
following form 

m n 

i=l i=l 

_ fo((log(l/e))'“+’”“e('“+’”“)"e('i+’"i)'*+'-l“l) li+mi>0 
|0((log(l/e))*“+’”“) Zi-|-mi = 0. 

Using the notation introduced above the proof is easy. We observe that each 
i G Lq and j G Mq contributes a factor log(l/e), whereas each i G Li, resp. 
j G Ml, gives rise to a factor resp. We define the tuples A', resp. 

Ii', by setting all components in A, resp. /i, which are less than s to 0. Then 
we obtain a total bound 0((log(l/e))*°“''”^“e^*i“''"‘i^'*“l^ I). If h + mi = 0 

then also |A'| -I- |^'| = 0 which proves the second case in (1^ . If -|- mi > 1 
we note that m — Z = |A| > |A'| -I- Iqs and |q;| — m = |p,| > |p,'| -|- mps. This 
implies — |A'| — |/r'| > (Zq + wo)s -I- Z — |a| and hence 1“!'^ I < £:(*o+iTio)s+i-|a|_ 
Inserting this into the above total bound matches the first case in (1^ and 
proves claim 2. 

We are now in a position to estimate the regularity r of u using (PI and 
(1^ . From PI we learn that r < t; and since s > 0 this is compatible with 
the assertion in la. In order to investigate the asymptotic behavior of (1^ if 
|a| > 1 we consider the cases s = 0 and s > 0 separately. 

s = 0: We recall that r = min(r, 1) — cr < 1 and we have to show that 
(I24II is Combination of (I25II (note that |a| > Z > 1) with 

the three possible cases 0(1), 0(log(l/e), 0{£*~’‘) of the growth rate of 
^1 directly yields an upper bound of the form 0(£™'"^*4)-|a|^jQg^2./£))^) 
(where k < Iq + mo -I- 1). Since the logarithmic factor is dominated by 
£~'^, for any cr > 0, the assertion is proved. 

s > 0: Now r = min(t, s-|-1) can be any real number and we have to go through 
all cases relating the possible values of |q!| and r. 

|a| < r: Since 1 < Z < |q;| the factor |6e^| is 0(1), and in (1^ . it^ we find 
|A| -I- 1^1 = |a| — Z < s, which in turn yields Zi = mi = Zq = mo = 0. 
Therefore we have an overall bound 0(1). 

|a| = r: If Z = |a;| then |A| -I- |^| = 0 and hence Zq -I- mo = 0 in ll^ which 
means 0(1) for this part. The factor \bi^^\ gives at most 0(log(l/£)). If 
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I < |a| then I < t and so is 0(1). Since |A| + |/.t| < s we deduce 
Iq + niQ < 1 and (ED) ensures an overall logarithmic bound. Hence we 
have an upper bound of logarithmic order in both (sub)subcases. 

|a| > r: We note that s > 0 implies (log(l/e))*‘>+"*°£l*‘>+"^“)'* = 0(1) what¬ 
ever the value of Iq + mo > 0. Therefore the first case in EH), li +mi > 1, 
always yields a bound 0(e®''‘*“l“l). 

I < t: The ^-dependent factor in 12411 is 0(1) and both cases in 12511 are 
dominated by 0(e'’+^“l“l). 

I = t: I contributes a logarithmic factor. The first case in I25II then gives 

0 (£s-i-t-|a|) Qj: which the part e® can be used to suppress this logarithmic 
factor; hence a bound is 0(e‘“l“l). The second case in (1^ yields an overall 
bound which is some power of log(l/e) and therefore clearly dominated 
by 

I > t: Here = 0(e*“*). Adding the factor according to the first case in 
12511 then gives a bound 0(e‘“l“l“''®) = 0(e*“l“l). On the other hand, using 
the second line in (1^ provides an overall bound 0{e* * • (log(l/e))*“+'"°). 
If I < |a| splitting off leaves an additional positive e-power to com¬ 

pensate for the logarithmic term. If Z = |a| we can again reason, like in 
earlier cases, that Iq + mo = 0. So, all branches of this (subsub)subcase 
lead to a bound 0(e‘“l“l). 

Collecting the results of all (sub)subcases we have established the asymp¬ 
totic upper bound 0(e™™(*'^+i)“l“l) of 12411 . 


□ 

The previous theorem indicates that we may expect a seismic wave to be 
about one degree smoother than the irregular medium variation if the source 
is prepared appropriately. In principle this would enable one to deduce from 
measurements of the wave an upper bound of the (global) medium regularity: 
first, estimate a strict upper bound of the wave’s Zygmund regularity r via 
wavelet analysis of the data; then the medium regularity cannot be better than 
r — I. 
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Appendix: Characterization of Zygmund regular¬ 
ity via continuous wavelet transform 

The proof to be presented below is a destination of methods and basic setups 
drawing from a variety of sources. We briefly sketch the basics of these as a 
preparation. 

(i) Zygmund classes can alternatively be defined by a discrete Littlewood-Paley 
decomposition (cf. m), also called dyadic resolution (e.g., in EH)- Let (fio = 

and for j G N put = J^j-i We have 

= (pj{^/2) and the support of ipj (j > 1 ) is contained in the annulus 
2^“^ < ICI < 2-^+^. By construction, the family {^j)j>o is a dyadic partition 
of unity: V’jiO — 1- Similarly, the equation V’jiD)u = u holds with 

convergence of the series in 

(ii) The classical Holder-Zygmund spaces can also be considered as the special 

cases in Triebel’s family of Besov-Hardy-Sobolev-type spaces (cf. ini 

Chapter 2, in particular 2. 6 .5/(1)]). These spaces are defined, for any s G M, by 
■ 6^,00 := jw G I := supj>o 2^11 (-0)^11 < oo}. The definition 

is independent of the particular choice of ip (cf. Ell 2.3.2]). 

(in) Both families of spaces, as well as C/, are realizations of the classical 

Holder-Zygmund spaces when s > 0. Therefore we clearly have B^ ^ = C/ 
in this case. In fact, equality holds for all real s: By EOl 2.3.8] (resp. [U 
Prop. 8 . 6 . 6 ]), for any r G M the operators (1 — A)*"/^ (resp. (1 — A)“’’/^) on 
map B^^^ (resp. C'/”’’) isomorphically into B^J^ (resp. C/); therefore we 
obtain B^^^ = C* Vs G M with equivalent norms and JuJ*. We refer 

to these spaces as Zygmund spaces of regularity s. In particular, we deduce that 
the definition of C/ is independent of the choice of (p. 

(iv) In Meyer’s book (cf. ^1 Chapter 3]) the Holder-Zygmund spaces are treated 
as special cases of Bony’s two-microlocal spaces (where s > 0, s' = 0, xq 

arbitrary). In fact, it is this point of view which is underlying the proof of the 
characterization via the (’continuous’) wavelet transform given in the following. 

Proof of Theorem^ Recall that c5^o(]R") is the subspace of J^(K.") consisting 
of functions with vanishing moments of all orders. Throughout the proof we will 
make use of the following fact which will allow us to balance vanishing moment 
conditions with regularity properties in occurring convolutions. 

Lemma 17. If / G .5^ with moments up to order m — I vanishing then one can 
find functions /a G (]«] = m) such that 

/= E 

\a\—m 


If, in addition, / G o5^o the functions fa can be chosen to be in 

This can be shown by adapting the proof of [HI Section 2.6, Lemma 12]. 
Concerning the notation of various constants in the estimates to follow we 
will use the generic letter C, with subscripts if we want to indicate dependence 
on certain parameters. 
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Part (i): Applying the above lemma to g we have J2\a\=m 9°'{{ga)e) 

and obtain 


Wgu{., e) =u*gs = (-1)'"£™ ^ (d^u) * {ga)e- 

\oL\—m 


Since G ™ with s — m < 0 and g^ G we have reduced the proof of 
0 to the task of estimating ||u * ge\\L°° where u G Cl with s < 0 and g G y . 

IfT > 1 arbitrary then u = (/3 (_D/T)u+/^ il}{D/T^udtjt (with J^'-convergence) 
and therefore we have for any £ > 0 fixed 

/ OO 

'iIj{D/ t)u * ge dt/t. (26) 

Let T > lie > T12 and estimate the two terms in (I26II separately. 

Since ip{D/T)u = ip{D)u + ijj{D/T)udt/t we deduce (recalling that we 
may assume s < 0) 


\\ip{D/T)u\\^^ < \\<f{D)u\\^, 


+ H{D/t)u\\^^dt/t 


< c(i 



< CP-’^ < 2-®C£L 


Therefore we obtain 

MD/T)u * g,\\^^ < MDlT)u\\^^ < Ce\ (27) 

To estimate the integrand in the second term of we assume t > T and 

choose Ip G T) with ip = 0 near 0 and '0 = 1 on supp('0). It follows that 
J-~^ip G oS^o and ip{D/t)u * ge = ip{D/t)u * ip{D/t)gg. 

Choose r G N such that r + s > 0 and apply Lemma El to obtain func¬ 
tions ipa, |a| = r, satisfying lF~^ipa G c5^o and tF~^ip = Y,\a\=r = 

E|a|=r •^■^(^“'0a)- Then ip{D/t)ge = E|a|=r /t){D^ g)e and since 

||0(Z3/t)M||^cxj < we have the estimate 


\\lp{D/t)u*ge\\L^ <t ’■£ ’■ XI IIV’(T’A)M*^a(T'A)(T>“5)e|lL- 

|Q;|=r 

<t-’’£“’' X \\'PiD/t)u\\^^\\ipa{D/t)(D°‘g)e\\Li 

\Qi\—r 

< Ct“(’'+'*)£“’' max \\'pa{D/t){D°‘g)^\\^i. 

\a.\—r 

We show that the appearing L^-norms have bounds independent of t and e. 

Writing ipa{D/t){D°‘g)g explicitly as a convolution and rescaling by t via 
substitution of the integration variable we have 

ipa{D/t){D°‘g)eix) = £“" J P~'^{ipa){y){D°^g){^ - dy. 

For any 1, the second factor in the integrand is bounded by C;(l-|-1| — -^1^)“* < 
2'-Ci{l + If P)"'(l -f If P)' < 2'C/(1 -f If P)“'(l -f \y\^y since te > X > 1. 
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Assuming I > n/2 and integrating also over x we finally obtain a bound for 
\\'4}a{D/t){D°‘g)e\\^i of the form C/,a f £“"(l + |f dx = Q^a f (l + lzl^)-‘ dz 
which is indeed independent of t and e. 

Taking the maximum of all bounds over |a| = r we arrive at the conclusion 
that for alH > T 

If i? > T arbitrary then 

pR pR 

I! / 'ipiD/t)u* gedt/t\\ < / \\ip{D/t)u* geW^^dt/t 

JT L°° Jt 

pR 

< Ce-'' / dt = 

Jt r + s 

When oo this upper bound tends to 7 ^(£T)“’' T“® < Ce®. This completes 
the proof of 0. 

Part (ii): 

Lemma 18. Let r > 0 and k G N such that k > r. Assume that hj {j G Nq) is 
a sequence of functions in C^(R”) with the property that there is i? > 0 such 
that for all fd G Ng with \f3\ < k 

\\dPh,\\^^<BV\^\. (28) 

Then the infinite series 

OO 

(29) 

3=0 

converges uniformly and defines an element in (R"). 

Proof. Since < B for all j the series is absolutely and uniformly con¬ 

vergent and defines a continuous bounded function h{x). Hence it is immediate 
that \\ip{D)h\\^^ < ||/i||ioo- It remains to estimate \\Pf3{D/t)h\\^^ 

for alH > I. We start by picking q G Nq such that 2'? < f < 2"?+^ and split the 
necessary summation according to 


\Pip{D/t)h{x)\ ^^P\-ip{D/t)hj{x)\ 

3=0 

g—1 OO 

= Y, 2-^^P\i3{D/t)hj{x)\ + Y 2-^"t^\i3{D/t)hjix)\ =: Si{x) + S 2 {x). 

3=0 3=q 

The terms in S 2 can be estimated as follows 

and hence S 2 {x) is dominated uniformly by a convergent geometric series. 
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To find a bound for S'i(a;) we apply Lemma ITTI and rewrite ip{D/t), as with 
i/j in the proof of part (i), in the form 'ip{D/t) = t~'^J2\a\=k ^a{D/t)D°'. Hence 

\ij{D/t)h,{x)\ < Y. 

\a.\—k 

and we obtain 


q-l q-l 

Si{x) < 

j=o j=o 


Oq(k-r) _ 1 

= r- < C. 

2k-r _ 2^ - 


Since t > 1 was arbitrary and the constants in the estimates are independent of 
q the lemma is proved. □ 


Lemma 19. If IT G Om(®.” x IR+) and satisfies with W substituted for 
WgU, then 

de 

Jo £ 

Proof. We show that the limit of := f^_N W{.,s) * ggdsfe, as N ^ oo, 
defines an element in C®. As used already in part (i) Tyemma, |l 71 imnlies g^ = 
£"* E|a|=m 9°‘iiga)e) and hence 


^ 5“ 

\(y\—m 


^ ^ e"* ^IT(.,e) * {ga)s'j 


.,(N) 


de. 


For any a G Nq with \a\ = m the map 9“: (7®+"* —> C® is continuous, hence it 
suffices to prove convergence of in (7®+™ (as N —> oo) for each such a. 

By a dyadic subdivision of the interval [2“^, 1] we find a corresponding series 
representation of in the form 


N-l 2-^ 

= Y £™"^IT(.,£) * {gc)e{x)de 

j=o 72-.J-1 

JV-i .1 

= V 2-^(-+®) / 2^®IT(.,2-^-,?) * igc)2-^r,ix)g^-^dg 

-^1/2 

N-l 

J=0 

where we have changed the variable e = 2~^g. Note that \\2Bw{.,2~^g)\\^^ 
< (777® < C independent of j. Therefore the sequence Vj^a satisfies the condition 
(I28II of Lemma Ci for any fc G N with k > m + s > 0 since 

\d'^v,,^{x)\ < [' |2^®W(.,2-^'77)*(a'^5„)2-.,(x)|2^'l'^l77’"-l'^l-id77 

7i/2 

<2^-|^l||2^-®W(.,2-^'77)||^„||a'^5„||^i r;™-l^l-id77 = C^,„2^l^l. 

Jl/2 
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Application of Lemma IIXI completes the proof. □ 

Lemma 20. Let W G Om(R" x R+) then 

de 

Ji £ 

Proof. Let R > 1 and put vr = f^W(.,e) * gede/e. Then vr is smooth, 
temperate, and converges weakly to some v G as i? —*■ oo (cf. Q). 

Clearly, any derivative d°'vR converges to d'^v then. But letting the deriva¬ 
tive fall on the factor inside the integral defining vr produces additional fac¬ 
tors When |a| is large enough to compensate for the polynomial growth 

of W{y,e) with respect to e this ensures absolute convergence of the classical 
integral. Hence for all |a| sufficiently large is smooth, yielding that v itself 
is smooth. □ 


To finish the proof of part (ii) we apply 0 together with m and obtain, 
with some polynomial p. 


N de , , de 

u= / WgU(.,£)*Pe -h / WgU{.,e)*ge -hp. 

Jo £ Ji £ 


The second term is smooth by Lemma 1*201 and the first term is of Zygmund 
regularity s by Lemma El It follows that u differs from an element in Cl only 
by some smooth function. 

This completes the proof of Theorem Q] 
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